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We suggest two realistic f(R) and one F(G) modified gravities which are consistent with local 
tests and cosmological bounds. The typical property of such theories is the presence of the effective 
cosmological constant epochs in such a way that early-time inflation and late-time cosmic accelera- 
£ — ■ tion are naturally unified within single model. It is shown that classical instability does not appear 

here and Newton law is respected. Some discussion of possible anti-gravity regime appearence and 
' related modification of the theory is done. 
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1. It is expected that modified gravity (for a review, see [l|) approach which may be related with string/M-theoryQ 
may help in the understanding of the dark energy origin. For instance, various modified f(R) gravities have been 
constructed and applied to the description of the late-time cosmic acceleration 0, 0, [H, [H, 0, H| with also the check 
of local tests. Even the form of modified f(R) gravity may be reconstructed from the known universe expansion 
history[9(. Hence, such gravitational alternative for dark energy may be considered as an alternative gravity theory 
. subject that it passes cosmological bounds and Solar System tests. 
O i' Recently, it has been proposed the class of modified gravities with the effective cosmological constant epoch 0, 
ITT1 . \l$L \l3L Il5l |. The very simple versions of such theories [HI [HI, [HI, EH with vanishing cosmological constant may 
describe successfully the ACDM epoch and pass the local tests/cosmological bounds. Nevertheless, the early universe 
epoch is not included there and deviations from Newton law at large scales may occur. 

In the present letter we suggest two realistic new models of modified f(R) gravity which may describe the early- 
time inflation and late-time acceleration in the unified manner, extending the earlier proposal of ref.Q. Moreover, 
such theories successfully pass the Solar System tests as well as simplest cosmological bounds and they are free of 
instabilities. Nevertheless, the first of suggested models may develop the anti-gravity regime, that is why its simple 
modification to avoid this problem is done. Finally, viable F(G) gravity which also unifies inflation with late-time 
• 1 acceleration and may pass the local tests/cosmological bounds is suggested. 

o 

S = -1 / d A x^—g{R + f{R)) . (I) 

• rH 

\g^F{R) - R^F'(R) - 9llv UF'{R) + V M V„F'(i?) = ~T {m)lxv . (2) 

Here F(R) = R + f{R) and T( m w is the matter energy-momentum tensor. 

Recently an interesting f(R) model has been proposed in [TTI | . In the model f{R) is given by 

m 2 ci (i?/m 2 )™ 

Jhs{R) = , p/ 2 ,n . , > (3) 

C2{R/m z ) +1 



2. The action of general f(R) gravity is given by 



K 

Here f{R) is an arbitrary function. The equation of motion in /(i?)-gravity with matter is given by 



which satisfies the condition 



lim fHs(R) = const , 

R— >oc 

hmf HS (R) = , (4) 
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The second condition means that there could be a flat spacetime solution (vanishing cosmological constant). The 
estimation of ref. ll| suggests that R/m 2 is not so small but rather large even in the present universe and R/m 2 <~ 41. 
Hence, 

t fn\ m 2 ci m 2 ci / R \ ~ n 

Jhs(R) 1 — \ — ' 5 

c 2 (4 \m l ) 

which gives an "effective" cosmological constant —m 2 ci/c2 and generates the late-time accelerating expansion. One 
can show that 

ff2 ^ ^ ( 70 km/s • pc) 2 ~ (lO- 33 eV) 2 . (6) 

C2 

At the intermediate epoch, where the matter density p is larger than the effective cosmological constant, 

. m 2 ci 

P > . (7) 

C2 

there appears the matter dominated phase and the universe expands with deceleration. Hence, above model describes 
the effective ACDM cosmology. 

Although the model [ll[ is very succesful (see, however, ref.[I3] for description of the Newton law at large scales 
here), the early time inflation is not included there. We suggest the modified gravity model to treat the inflation and 
the late-time accelerating expansion in a unified way. We now consider simple extention of the model to include 
the inflation at the early universe. In order to generate the inflation, one may require 

lim f{R) = -A, . (8) 

R — >OG 

Here Aj is an effective cosmological constant at the early universe and therefore we assume Aj > (l(T 33 eV) . We 
may assume Aj ~ 10 20 ~ 38 . Instead of (JSJ), as in Starobinsky's model, one may require 

lim f(R) cx R 2 . (9) 

R — >oo 

In order that the accelerating expansion in the present universe could be generated, let us consider that f(R) could 
be a small constant at present universe, that is, 

f(R ) = -2R , f'(R ) ~ . (10) 

Here Rq is current curvature Rq ~ (lO~ 33 eV) 2 . The next condition corresponding to the second one in ([J| is: 

hm/(i?) = 0. (11) 

In the above class of models, the universe starts from the inflation driven by the effective cosmological constant © 
at the early stage, where curvature is very large. As curvature becomes smaller, the effective cosmological constant 
also becomes smaller. After that the radiation/matter dominates. When the density of the radiation and the matter 
becomes small and the curvature goes to the value Rq pI7|) . there appears the small effective cosmological constant 



(|10p . Hence, the current cosmic expansion could start. 
An example satisfying ([5]), (fTU|) . and (TTTj) is 

f(R) = (R-Ro) 2n+1 + R 2 o n+1 — ^ = i VA (12) 

f + h {(R Rof n+1 + Rl n+1 } ^ fo + h{(R- Ro) 2n+1 + Rl n+1 } 

Here n is a positive integer, n — 1, 2, 3, • • • and 

R 2 n+1 ^ 1 



h + hRl n+1 /: 

that is 



2R , - = /a , (13) 
i 



p2n p2n+l n2n -i 

f °-~2 A^ 2~ ' fl -A~- (W) 
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By introducing the auxilliary field A one may rewrite the action ((T|) in the following form: 

S = \ [ d 4 xV^{(l + f{A)) (R-A)+A + f(A)} . (15) 



From the equation of motion with respect to A, it follows A = R. By using the scale transformation g^ v — > e° ' g^ v 
with a = — In (1 + f'(A)), we obtain the Einstein frame action: 



A F(A) 



(16) 



S ° = ^ j dix ^\ R -l{W)) ^"^-Aa) ' F W 
= J d^x^-g (r - Z -gP"d p ud a a - V(a)) , 

V(a) = e°g ( e -) - ff (g (e"')) = ^ - • (17) 

Here g(e~ a ) is given by solving cr = — In (1 + f'(A)) = InF'(A) as A = g(e~ a ). After the scale transformation 
gp.ii — * e° r fl l J ui/, there appears a coupling of the scalar field cr with the matter. For example, if the matter is the scalar 
field $ with mass M , whose action is given by 
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\ J d*aV=ff - M 2 $ 2 ) , (18) 

there appears a coupling with cr in the Einstein frame: 

S*e = \J dtxy/^g^ffcrdy&du* - M 2 c 2ct $ 2 ) . (19) 
The strength of the coupling is of the gravitational coupling k order. Unless the mass of cr, which is defined by 

a _ 1 dPV(cj) _ 1 J A AF(A) 1 1 

CT_ 2 da 2 2|F'(A) (F'(A)) 2 J 1 j 

is large, there appears the large correction to the Newton law. 

In air on the earth, the scalar curvature could be given by A = R ~ 10~ 50 cV 2 . On the other hand, in the solar 
system, we find A = R ~ lCT 61 eV 2 . In the model (JUJl with (13]), if n > 3 in the air or n > 10 ~ 12 (if A, ~ 10 20 ~ 38 ) 
in the solar system, we find 

/o « h{(R- Ro) 2n+1 + R 2 n+1 } ~ hR 2n+1 (21) 



and 



F{R) = R+m~R-y i + 1 t^~j x - (22) 

Furthermore, if n > 6 ~ 7 (if A, ~ 10 20 ~ 38 ) in the air, we find F'(i?) - 1 and 

™^-||^-^ 10 38 ^ 6 eV 2 , (23) 

which is very large and there is no observable correction to the Newton law. We should note 1 mm ~ (10 _4 eV) _1 . 
On the other hand, in the Solar System, if n > 10 ~ 12 (if A, ~ io 20 ~ 38 ) , 

m 2 _ I ^° _ 1Q 239~295-10n y2 (2A] 

° 2F"(A) hR^+i U 6V ' 1 ' 

which is large enough since the radius of earth is about 10 7 m ~ 10 _14 eV) and there is no visible correction to the 
Newton law. 
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There may exist another type of instability (so-called matter instability) in f(R) gravity [T(|(see, also [13. fl4|). 
The instability might occur when the curvature is rather large, as on the planet, compared with the average curvature 

at the universe R ~ (l0~ 33 eV) . By multipling Eq.([2|) with g^ v ', one obtains 

np F^{K)_ pv7pp F%R)j 2F(i?) K 2 

+ FW(R) p 3FW(R) 3F(2)(i?) ~ 6F( 2 )(i?) " w 

Here T = T,£* and F^ n \R) = d n F(R)/dR n . We consider a perturbation from the following solution of the Einstein 
gravity: 

R = R b = --T > . (26) 

Note that T is negative since \p\ <C p on the earth and T = —p + 3p ~ —p. Then we assume 

R = Rb + R P , (\Rp\ « \R b \) • (27) 

Now one can get 



FW{R b ) " 3F( 2 )(i? b ) 3FW(Rb) 3FW{R b ) 
F^(Rh) 

+DR P + 2 f( 2)\^ p r ^ pr p + U(R b )R P , 
( F^{R b ) F(*){R b ? \ R b 

{WkrT)- FW(R b )*) vpRbV Rb + Y 

F^(R b )F^(R b )R b F^(R b ) , 2F(R b )F^(R b ) F^(R b )R b 



(28) 



3F( 2 )(J? 6 ) 2 3F( 2 )(i? h ) 3F( 2 ) (i? b ) 2 3F( 2 ) (i? b ) 2 ' 

Since Di? p ~ —d 2 R p , Eg. ([28)1 has the following form: 

= -<9 2 i? p + U{R b )R p + const. . (29) 

Then if U(R b ) is positive, R p becomes exponentially large as a function of t: R p ~ qs/u^) 1 an( j th e system becomes 
unstable. In the model (fT2"|) with (p~5|) . if n > 2, we hnd 

, , (2n + Z)fiRl n+2 , , 

3(2n + l)/o 

Therefore there is no this kind of the instability. 

Hence, in the model (fT2")l . the universe could start from the inflation. As curvature becomes smaller, the effective 
cosmological constant becomes small and after that the radiation/matter could dominate. When the density of the 
radiation and, later, of the matter becomes small and the curvature goes to the value Rq, there could appear the 
small effective cosmological constant and the late-time accelerating expansion starts. 

3. As is clear from (|15l) . if F'(R) = 1 + f'(R) > 0, the square of the effective gravitational coupling becomes negative 
K 2 ff = k 2 /F'(A) and theory enters the anti-gravity regime. Nevertheless, if the theory enters such regime at the 
future, such theory may still be viable. We now check if the model (|12p could pass the anti-gravity constraint. Since 



f'(R) 



(2n + l)/o (-R - Rof n 



(h + h{(R-Ro) 2n+1 + R^ +1 }y 



2(2n + l)/ (R - R a ) 2n - 1 ( -/„ - + h(n + 1) (R - R Y 
f"(R) = * 3 '- , (31) 

(fo + fi{(R-Rof n+1 + R 2 n+1 }y 

f'{R) has a maximum at R = i?, which satisfies 

(R-Ro) = /0 f + / " -tAtT- (32) 
V / 1 1 In + 1 1 1 (n + 1 
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Then the maximum value of f'(R) is given by 



2n+l 



f'(R) ~ ~^r- ~ (Rofx)-^ « -1 • (33) 

Here we have used (fHj) . Therefore F'(R) < and the anti-gravity regime could occur in general. We can consider 
only the region where F'(R) > or f'(R) > —1 without the the anti-gravity problem. From (|32|) , however, one gets 

R~ (lO" 33+ ^ev) 2 , (34) 

which could be small compared with the curvature at the early universe. There are two solutions for the equation 
f'(R) = — 1. Let denote the larger one by R+ and smaller one by R-. By assuming 

f « h { (R+ - Ro) 2n+1 + Rl n+l } ~ hR 2n+1 , (35) 



we find 

R, 

On the other hand, by assuming 
we find 



(2n+l)V 1/(2 " +1) 



,.- , - (lO^^+^r^ev) 2 . (36) 

fx J ^ ' 

fo » A {(R- - R ) 2n+1 + R% n+1 } , (37) 



R_ ~ R ~ (lCT 33 eV) . (38) 

Then R± could be rather small compared with the scale of the inflation (l0 19 ~ 25 eV) 2 . Then the model (p~2|) itself 
seems to be not viable. It could be possible to replace f'(R) given from (fl2|) in a region R_ — e_ < R < R+ + e+, 
where e± are small positive constants, with a proper function whose first derivative is always greater than — 1. For 
example, 

f /oid(-R) R < R- - e- 

f(R) = I f id(R--e-)+fou( R -- e -)( R - R -+£-) R--e-<R<R+e+ . 

{ /oid(-R) + /oi d (#- - e-) + e+-R-+ e_) - /oid(i?+ - e+) + /old(#- - e_) i? > 

/ p p \2n+l . p2n+l 

/om - — {R r o) + ° = • (39) 



/o + / 1 {(i?-i?o) 2n+1 +i? 2 " +1 } 



We now choose e± to be /^ ld (i?_ — e_) = f^ ld (R + + e+). Hence, one obtains /'(-R) > —1 and there is no problem of 
anti-gravity so such corrected model seems to be quite realistic. 

From another side, in order that the model (fT3>|) could be viable, there should occur the discrete transition from 
R = R + and R = R- . Such a transition might occur as in the usual phase transition with the jump of the value of 
the order parameter. The properties of such possible transition should be studied in detail. 

4. In order to avoid the above anti-gravity problem from the very beginning, we may propose another model given by 



f{R) = -fo / dRe c*-"!)*" " A - . (40) 
Jo 

Here a, (3, fo, and R\ are constants. Then by construction, as long as < /o < 1, f'(R) > — 1 and therefore 
F'(R) > 0. Since 



pRi oaf" 
f(Ri) ~ -fo / dRe Tn^no™ = -f Q A n (a)R 1 , 
Jo 



A n (a) = / dxe ^ , (41) 
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and — f(Ri) could be identified with the effective cosmological constant 2i?o, we find 

foAni^Ri = R . (42) 
Note that A n (0) = 1, A n (+oo) — 0, and A'(x) < 0. On the other hand, since 

/>oo 

/(+oo) ~ / dRe- R ^ = -/oMi , (43) 



and — /(+oo) could be identified with the effective cosmological constant at the inflationary epoch, Aj, we find 

faP = 1 • (44) 
Let us now investigate the correction to the Newton law. One gets 

aR i n r / On E?2n i \ 

f"(R) = -f n e" — --i — ) . (45) 

If in air on the earth n > 3 ~ 4 (if Aj ~ 20 20 ~ 38 ) or in Solar System n > 10 ~ 12, we find m 2 (TSTTf) is given by 



in 



2 - ' 0A * 'l0 10 ~ 16 GeV) 2 , (46) 



2F"(R) 2/ 



which is very large and positive and therefore the correction to the Newton law is very small. 

We now also investigate the instability indicated in [10]. Inside the earth, U(R) (|28[) has the following form 

U{Rb) ~ ~ "X < ° ■ (4?) 

Here we have assumed n > 10 ~ 12. Since U(Rb) is negative, there does not occur such the instability. 

Hence, in the model J30|), the universe could start from the inflation driven by the effective cosmological con- 
stant — /(+oo) at the early stage, where curvature is large. As curvature becomes small, the effective cosmological 
constant becomes small too and radiation/matter dominates at the intermediate epoch. When the density of the ra- 
diation/matter becomes small and the curvature goes to the value R\, there appears the small effective cosmological 
constant — /(i?i), so that the accelerating expansion starts. Hence, the realistic modified gravity consistent with Solar 
System tests and unifying inflation with cosmic acceleration is constructed. 

Instead of /(i?)-gravity, we may consider the F(G)-gravity, where the action is given by [l6| 

S=\ ( d 4 x^j(R + F(G)) , (48) 



<- 2 

t\ J 

and G is the Gauss-Bonnet invariant: 

G = R 2 - AR^R^ + R^R^P" . (49) 

Note that in F(G)-gravity, there are no problems 16] with the Newton law, instabilities and the anti-gravity regime. 
One may consider the model similar to (JT2J) : 

P-W^C i Mi T . ,50, 

Here Go corresponds to the present value of the Gauss-Bonnet invariant. Since F'(G) = when G = Go and G = +oo, 
F(G) becomes almost constant and can be regarded as the effective cosmological constant. As in ([15)) , we may identify 
F(oo) as the cosmological constant for the inflationary epoch and F(Rq) as that at the present accelerating era. Then 
instead of HU), one gets 

F °-~2 A, 2~ ' Fl -Tr (51) 



Hence, the universe starts from the inflation driven by the effective cosmological constant Aj (f5"Tj) at the early epoch, 
where curvature and therefore the Gauss-Bonnet invariant G is very large (for the study of cosmological perturbations 
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in such theory, see [I3|)- As curvature becomes small, the effective cosmological constant becomes small too and 
radiation/matter could dominate at the intermediate universe. When the density of the radiation/matter becomes 
small and G goes to the value Go hi |5D|) . there appears the small effective cosmological constant —F(Gq), and cosmic 
acceleration starts. 

Thus, we suggested several realistic models of f(R) and F(G) gravity which propose natural unification of the 
early-time inflation and late-time acceleration being consistent with local tests and cosmological bounds. Definitely, 
more precise local tests/cosmological checks for such theories should be made in future, having in mind that more 
precise observational data will be available soon. 
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Appendix We now consider how the universe can reach the exit of the inflation. For simplicity, we consider the 
model in (|12p . In order to compare with usual scenario of the inflation, we work in the scalar-tensor form in the 
Einstein frame (fit?)) . In the epoch of the inflation, the curvature R — A could be large and we find 

/(*)~-7- + 7Ji&TT. ( 52 ) 



h flR 2n+1 



and therefore 



2n + l 



2n + l)/ v{a) „ 1 _ 2(n+l)/ / fj*\ ^ ^ 



A 2 A 2 "+ 2 ' w h ft V(2n + l)/o 

Since a is now dimensionless, the condition for the slow roll could be given by | V /V\ <C 0. Now we have 

V>) - dJ^.Y*'. (54) 



(55) 



V(a) Jl \(2n+l)f 
If we start with a ~ 1, by using (|14l) . we find 

V'(a) (Rl n 



V(a) \A 



t 2/i 



which could be very small and the slow roll condition could be satisfied. The potential V(a) in (fS"3"|) tells that if we 
start with a ~ 1, the value of a increase very slowly and therefore R becomes smaller. Then a becomes large enough, 
the inflation could stop. 
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